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Abstract 

We use time-dependent density functional theory to calculate the energy loss of an 

antiproton colliding with a small Al cluster previously excited. The velocity of the 

antiproton is such that non-linear effects in the electronic response of the Al cluster 

are relevant. We obtain that an antiproton penetrating an excited cluster transfers less 

energy to the cluster than an antiproton penetrating a ground state cluster. We quantify 

this difference and analyze it in terms of the cluster excitation spectrum. 
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1. Introduction 

A charged particle moving across a metallic target is able to create electronic ex- 
citations in the medium at the expense of its own kinetic energy. Research on this 
phenomenon has been broad in condensed matter physics and materials science because 
of its relevance in various fundamental and applied topics, such as radiation damage, 
medical physics, and ion sputtering. 

A key point in the theoretical analysis of the slowing down of charged particles in 
metals is the intensity of the perturbation that the moving particle introduces in the 
medium. For a particle of charge Q moving with a velocity t; in a standard metal, 
the perturbation strength can be roughly characterized by the Sommerfeld parameter 
Tj = Q/v [IJ. If such ratio is small, rj << 1, linear theory is naturally applied and 
accurate results for the particle energy loss are found. U rj » 1 and v is much smaller 
than the typical velocities of the electrons in the medium, various non-perturbative 
methodologies have been successfully applied PJ. In between these two cases, in the 
regime of intermediate velocities, accurate descriptions of the energy loss process are 
much more involved because quasistatic or perturbative approximations break down 
even for unit-charge projectiles. Only recently calculations based on time-dependent 
density functional theory (TDDFT) PHEIEIITIEIEIIIO] have shown its potential to 
close this gap. 

The achievements of TDDFT in the non-linear description of electronic excitations 
pave the way to answer new and challenging questions in the field. In the traditional 
description of energy loss processes, the target is always considered as initially in its 
ground state. However, the energy lost by a travelling charge in a metallic medium 
should be affected by the electronic state in which the target is. Otherwise said, if 
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electronic excitations have been already created in the system, the electronic response 
to the incident perturbation, and consequently the energy loss, will be different. In this 
work we try to quantify this difference for the particular case of a point charge crossing 
metallic clusters of a few A size. 

The experimental investigation of the excitation and ionization of neutral metal 
clusters by collision with positively or negatively charged particles has been intensive. 
In particular, the ionization of metal clusters by low energy singly and multiply charged 
ions and protons [II1|12] and the ionization of neutral metal clusters by slow electrons [I3l 
[TU [ISl US] have been studied. Description of such processes from the theoretical point of 
view is incomplete and requires further investigation. In our work we study the collision 
of an Al cluster with a slow negative point charge (an antiproton). The choice of the 
antiproton as a projectile allows us to avoid complications related to the electron capture 
by the cluster if the projectile is an electron and the electron capture by the projectile if 
the latter is an ion. Our goal is to identify the distinct effects that arise in the dynamic 
screening and the projectile energy loss when the metallic target has been previously 
excited by a preceding projectile. In spite of the fact that our model is simplified, the 
results of our study can contribute to the understanding of the fundamentals of the 
dynamical processes during collision of charged particles with metallic clusters. 

We perform an explicit time propagation of the electronic state of the system using 
TDDFT and evaluate the energy lost by the charge when crossing ground-state clusters. 
We compare this quantity with the amount of energy lost when the projectile crosses a 
cluster excited from a previous collision. We show that the difference is appreciable and 
give the explanation of this change as a consequence of the excited state of the cluster 
as well as of the emission of electronic charge from the excited cluster. 

Non-linear effects in the excitation of metal clusters have been previously analyzed 
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with TDDFT. In particular, electron dynamics in clusters under intense laser fields 
are an active hot topic of research [T7j because of the possibilities offered to explore 
and control ultrafast processes. The resonance energy of collective excitations in these 
systems has been shown to depend on the intensity of the perturbation |T8]. Here we 
focus on a different type of external perturbation, namely, that derived from a point 
Coulomb charge crossing the system. We will show, nevertheless, that similar shifts in 
the position of the plasmon peaks are found. 

Hartree atomic units {e = h = rrie = I) will be used throughout this work unless 
otherwise stated. 

2. Methodology 

Let us first define the system under study: We will focus on a negative point charge 
(an antiproton) crossing a metal cluster. Electron dynamics in metallic systems typically 
lie in the femtosecond and subfemtosecond time scales [121 [2D]- For this reason and for 
the kind of processes that we study, we consider the cluster ion cores as frozen. We 
further simplify the problem and use the spherical jellium model (JM) to represent the 
cluster. In the JM, the ions are substituted by an homogeneous background of positive 
charge with density t^q (r) = nQ{rs)Q{Rci —r). Here Rd is the radius of the cluster, Q{x) 
is the Heaviside step-function and nQ{rs) is the constant bulk density, which depends 
only on the Wigner-Seitz radius (1/?t-o = 47rr^/3) [21] . The number of electrons in a 
neutral cluster is = {Rci/rsf ■ 

The ground state electronic density of the cluster n{j) is obtained using the Kohn- 
Sham (KS) scheme [22] of density functional theory (DFT) [23j. The ground state KS 
wave functions are expanded in the spherical harmonics basis set [21]. 
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The evolution of the electronic density in the cluster in response to the field of the 
moving charge is calculated using TDDFT [25]. We propagate the ground state wave 
functions ipi{r,0) = (pi{r) solving time-dependent KS equations: 



The effective potential includes four terms Veff{r,t) = Vea.t(r,t) + VH-(r,)f:) + Vxc{r,t) + 
Vp{r,t), where Vext is the external potential created by the positive background. Vh is 
the Hartree potential created by the electronic density. Vxc is the exchange-correlation 
potential, calculated with the standard adiabatic local density approximation (ALDA) 
with the Perdew-Zunger parametrization of Ceperley-Alder exchange and correlation 

potential [26]. Finally, Vp{r,t) = 9(t) is the potential created by 

the antiproton and acting on the valence electrons of the cluster. We use cylindrical 
coordinates (p, z) in the time-dependent calculations, which are more appropriate since 
the problem has axial symmetry. The origin of coordinates is located at the center of the 
cluster. The antiproton is represented by a negative point charge {Qp = — 1) which moves 
with constant velocity v along the 2;-axis. At time t=0, the antiproton is located at a 
distance from the cluster (50 a.u.) far enough to avoid a significant interaction between 
the projectile and the target. The time propagation of the electron wave function is 
performed using the time-stepping algorithm: ipi{p,z,t + dt) = e~^^^^^ipi{p, z,t). The 
split operator approximation is then used to separate the potential and kinetic energy 
terms in the e"*'^*'^* time propagator. The action of the kinetic energy operator is 
calculated using Crank-Nicolson propagation scheme. A detailed description of the 
numerical procedure can be found in Refs. [27t |28| [29] . 

From the time-dependent KS orbitals we obtain the time-evolving electronic density 
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of the excited cluster n{p,z,t) — \ipi{p, Zjt)]"^. The force acting on the moving 

iSocc 

antiproton along the z-axis is obtained from the time-dependent electronic density and 
includes the effect of the positive background: 

^ 2./, dp a.^^^i0^0^^[Mt) - 4 (2) 

The energy loss is then obtained from the integral: 

POD 

Eloss ^-v F,{t) dt. (3) 
Jo 

We will study the energy loss in two different motion cycles. In the first cycle, the 
antiproton moves towards the cluster with a constant velocity v, crosses it following a 
symmetry axis through the cluster center, and eventually moves away until it reaches a 
turning point arbitrarily defined. The turning point is at a distance from the cluster far 
enough not to have any residual interaction. The cluster is then left in an excited state. 
The electronic energy transferred to the cluster during the collision is calculated. In the 
second cycle, the antiproton turns back from the turning point and starts to approach 
the excited cluster with the same constant velocity v. In the second crossing of the 
cluster, the latter now in an excited state, energy is again transferred to the cluster. We 
calculate the energy lost in this second cycle and compare the obtained value with that 
of the first cycle. 

3. Results and discussion 

We have chosen a small Al (r^ = 2.07) cluster with N =18 electrons and with radius 
Rci — 5.43 a.u. (~ 0.29 nm). In all the calculations shown in this article, the projectile 
velocity is v — 0.5 a.u. The ALDA-TDDFT method used here predicts very well the 
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energy loss of antiprotons in Al targets. The method gives very good agreement with 
measurements in Al bulk for antiproton velocities up to 1.8 a.u. Above this velocity the 
excitations of the inner shells in Al start to contribute to the energy loss and results 
deviate from the experimental ones 

The antiproton starts its motion at time t = from the position zq = —50 a.u. 
After the first collision, the projectile continues to move until t = 1000 a.u. At this 
time the second cycle starts and the antiproton takes the way back to collide again 
with the cluster. We call r to the time interval between both collisions. In both cycles 
we calculate the force Fz{t) experienced by the projectile due to the interaction with 
the cluster through Eq. [2] From Fz{t) we obtain the value of the energy lost by the 
antiproton Eioss by means of Eq. [3j In addition, we consider three other different time 
spots for the second cycle to start: 1003.5, 1005, and 1010 a.u. The purpose of using 
different time delays is to check the sensitivity of the final result to the dynamics followed 
by the electron density in the cluster excited state. With our choice of time delays, the 
antiproton reaches the excited cluster respectively Ar = 7, 10 and 20 a.u. of time later 
than in the reference calculation. Depending on the value of Ar, the antiproton will 
start to cross the surface of the excited cluster meeting a minimum or a maximum in 
the electronic density oscillations, or an intermediate state. The density oscillations will 
be discussed later. The results for the energy loss are summarized in Table |T] 

The first interesting conclusion that can be extracted from the results of Table [T] 
is that the energy loss of the antiproton crossing the excited cluster is consistently 
lower than the corresponding value for the antiproton colliding with the cluster in the 
ground state. There are two reasons for the decrease of the energy loss. One reason is 
that, after the first collision, the cluster is emitting an amount of electronic charge that 
roughly corresponds to one electron. This means that during the second collision the 
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antiproton is interacting with a smaller amount of electronic charge and therefore loses 
less energy. In order to check the relevance of the change in the electronic charge of 
the cluster, we performed an additional calculation, namely that of the energy loss in a 
positively charged cluster which contains 17 electrons and remains in its ground state. 
The obtained value of 0.8211 a.u. is lower than the value of the energy loss in the neutral 
cluster, which is given in the Table |T] and is equal to 0.8527 a.u. The difference between 
these two results is around 4%. However, the difference is not as big as in the case of 
time delays Ar = and 20 a.u. given in the table, which is up to 11% of the value of 
the energy loss in the first collision. This allows us to conclude that the emission of one 
electron from the cluster only partially explains the observed decrease of energy loss in 
the second collision. 

Another reason for the lowering of the energy loss is that the cluster is excited after 
the first collision with the antiproton. Namely, the electronic density of the cluster is 
starting to oscillate in time with a given frequency. As we mentioned before, depending 
on the value of the time delay Ar, in the second collision the antiproton meets different 
states of the electronic density oscillations at the surface of the cluster. In what follows 
we are going to analyze the difference in the energy loss between two collisions depending 
on the time delay of the second collision. 

For different time delays Ar of the second collision, the value of the energy loss 
slightly varies. In order to illustrate this, we show the difference in the force between 
the first and the second collision AF, = F^'^* — .p^nd £ different values of Ar. Here 

^ Z Z 

is the force felt by the moving charge colliding with the non-excited metal cluster, 
which is equal for all r; _p2nd force felt by the antiproton colliding with the excited 

cluster. AFj. is shown in Fig. [T] as a function of the antiproton position. In this figure, 
large negative values of indicate positions of the antiproton before each collision with 
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the cluster. The total force during the first collision F^^^ is shown in the inset of Fig. [T] 
as a function of the projectile position. The two strong features in F^'"^^ correspond to 
the antiproton crossing the cluster surface. Away from the cluster, the antiproton is 
attracted by the induced dipole. Inside the cluster, the electronic density rearranges in 
order to screen the strong perturbation created by the moving antiproton. The force 
inside the cluster oscillates about a mean value that roughly corresponds to the effective 
stopping power for this particular velocity of the projectile {v = 0.5 a.u.) The curves 
in the main panel of Fig. [T] show how the force felt by the antiproton changes depending 
on the time at which the second collision starts. 

The fact that the energy loss is different at different time delays of the second collision 
can be also analyzed looking at the total energy of the cluster. The energy is shown in 
Fig. [2] as a function of the antiproton position. From Fig. [2] we can see that the total 
energy of the cluster is increased by the collision. This increase in energy is the value 
of the energy transferred by the antiproton to the cluster or, in other words, the energy 
lost by the antiproton. We can see as well that, in all cases, the energy loss after the 
second collision is lower than after the first collision. The curves for At =0 and for 
At =20 a.u. are similar. This is consistent with the values of the energy loss given in 
Table [T] for these two cases. We can also see the longer range of the cluster-antiproton 
interaction during the second collision. This is due to the net positive charge of the 
excited cluster. 

The dependence of the energy loss on the time delay between collisions can be under- 



stood by looking at the time evolution of the induced electronic density. Figure 3a shows 
the change in electronic density An{z,p = 0,t) = [n{z,p = 0,t) — n{z,p = 0,t = 0)] 
inside the cluster in units of the background density uq, along the z— axis and as a func- 
tion of time. The results are shown for the calculation with At =0. The time interval in 
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Fig. 3a is chosen to include the moment at which the second colhsion of the antiproton 
with the cluster takes place. In the figure, the projectile moves from the right to the 
left. The second collision starts at t ^ 1890 a.u. We clearly see the negative change in 
density, originated by the Coulomb repulsion between the incident antiproton and the 



cluster electrons. From Fig. |3a] we can see that the excitation created by the moving 
charge in the cluster leads to oscillations in the induced electronic density: Minima and 
maxima in the induced density are observed. Depending on the time delay between 
collisions r, the impact of the incoming antiproton with the previously excited cluster 
can bump into a minimum or a maximum of the electronic density oscillations. In the 
first calculation (Ar = 0) and when the time delay is Ar = 20, the antiproton starts to 
cross the excited cluster when there is a maximum in the electronic density oscillations 



at the surface of the cluster (the change in density in Fig. 3a is positive). In the case of 
Ar = 10 a.u. the second crossing finds a minimum of the electronic density oscillations 
at the cluster surface. The time delay Ar = 7 a.u. is chosen to have a case in which the 
second crossing falls neither on the minimum nor on the maximum of the change of the 
electronic density, but in-between these two situations. Depending on this circumstance, 
the value of the energy lost by the antiproton varies. 



From Fig. 3a, we can also see that the minima and maxima in the induced electronic 
density become more pronounced after the second collision, indicating that the cluster 



is further excited by the second collision. This can also be seen in Fig. [3b] where the 
change in density is shown as a function of time for the particular value of z = 4 



a.u., marked with a dashed line in Fig. 3a The amplitude of the electronic density 
oscillations increases after the second collision. This is also observed in Fig. |3c] and 



Fig. 3d where we illustrate the density distribution in the cluster before {t = 1812 a.u.) 
and after (t = 1992 a.u.) the second collision. The change in density is plotted in a 
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plane in (p, z) coordinates with the center of the cluster at {p = 0,z =0). The negative 



and positive peaks in the right panel (Fig. 3d) are much more intense than in the left 



panel (Fig. 3c). These pronounced oscillations show that, after the second collision, the 
oscillations of the induced electronic density are stronger. The excitation created by the 
second antiproton enhances that created during the first collision. However, the similar 
distribution of the induced charge seems to indicate that similar electronic modes are 
excited in both events. 

In order to calculate the frequency of the density oscillations we perform a Fourier 
analysis of the time evolution of the dipole moment P{t) — )■ P{uj) created by the elec- 
tronic density in the excited cluster. The Fourier transform is done for two different 
cases! a ) after the single collision and without including the second collision, and b) 
after the two collisions. In this analysis we use the time evolution of the dipole during 
~1200 a.u after each collision and an exponential mask function (centered in the middle 
of such interval) to avoid spurious effects due to the use of finite time interval. The 
results for the respective dipole power spectra |P(a;)p are shown in Fig. |4j Two peaks 
are shown by arrows at frequencies u = 0.261 and 0.284 a.u. (corresponding periods of 
the plasmon oscillations are T ^ 24.1 and 22.1 a.u.). They roughly correspond to the 
expected value of the plasmon energy in the cluster: The plasmon frequency for a per- 



fect metal sphere can be calculated from the value of the density hq as Up = a/47™o/3 
(Mie plasmon frequency) [30] . For an Al cluster (r^ = 2.07) this value is Up = 0.34 
a.u. In our calculations the obtained frequency is lower than the frequency given by 
the classical Mie theory. This is due to the small size of the cluster and because we 
use quantum theory for the calculation of the frequency. A red shift with respect to its 
classical Mie value is frequently observed in clusters of simple metals [311 ESI El] • 
In Fig. |4| we also see that the plasmon peak is shifted to higher frequencies after the 



11 



second collision. This behavior is consistent with the non-linear shift of the plasmon 
frequency under non-perturbative conditions [18]. The blue shift of the frequency af- 
ter the second collision is in part related to the emission of electronic charge from the 
cluster due to the interaction with the antiprotons. It was observed that the resonance 
position moves to higher frequencies when increasing the positive charge of the cluster 
|32j . Lower-energy excitations are also present in the spectra and can be attributed to 
excitations of electron- hole pair character [T8] . 

4. Conclusions 

In summary, we have calculated the energy loss of an antiproton colliding with a 
small Al cluster, both when the cluster is in the ground state and when the cluster is in 
an excited electronic state. We have shown that the antiproton loses less energy when 
penetrating a cluster previously excited. The lowering of the energy loss is related not 
only to the fact that the cluster is transferred to an excited state, but also to the fact 
that the cluster loses one electron during the first collision with the antiproton. 

We have also shown that the projectile creates a plasmon in the cluster and that the 
plasmon peak shifts to higher frequencies in the second collision. This corresponds to 
the observed shorter period and larger amplitude of the electron density oscillations in 
the cluster after the second collision of the antiproton with the cluster. The shift of the 
plasmon peak to higher frequencies is partially due to the emission of one electron from 
the cluster, which thus becomes positively charged. 

Our work is another example of how TDDFT in the time domain is an extremely 
useful tool to study electron dynamics in finite-size objects, as well as to analyze the 
energy loss processes of charges interacting with condensed matter. 
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Table I: Energy loss Eioss (in a.u.) of an antiproton crossing the spherical Al cluster in 
ground and excited states. 





1^* collision 


2nd coUision 


2nd collision, 

At = 7 au 


2nd collision, 

At = 10 au 


2nd collision. 
At = 20 au 


^loss 


0.8527 


0.7583 


0.8318 


0.8099 


0.7554 
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Figure 1. Difference in the force between the first and second coUision of the antipro- 
ton with the cluster, AF^, for different time delays r between collisions as a function of 
the projectile position Zap. Inset: total force during the first collision F^**, as a function 
of the projectile position. Dashed hues show the borders of the cluster (Rd — 5.43 a.u.). 
All quantities are in a.u. 

Figure 2. Total energy of the cluster Efot for different coUisions and for different 
time delays At, as a function of the antiproton position Zap- AH the energy curves 
corresponding to the second crossing are referred to the value of the energy prior to the 
first crossing when antiproton is far from the cluster. Dashed lines show the borders of 
the cluster {Rd — 5.43 a.u.). All quantities are in a.u. 

Figure 3. (a) Time evolution of the induced electronic density inside the cluster along 
the axis (p = 0.02 a.u.) including the time at which the antiproton crosses the excited 
cluster. The color code shows the change in density [n{z, p = 0,t) —n{z, p = 0, t = 0)] in 
units of the background density Uq. The dashed line in panel (a) indicates the position 
z — A a.u. for which, in panel (b), we show the change in density as a function of time. 
Dashed line in panel (b) indicates the moment when the second collision starts, (c) and 
(d) show the change in the electronic density [n(r, t) — n{r, 0)] of the spherical cluster 
(color codes) in (p, z) coordinates at times t — 1812 a.u. and t — 1992 a.u. respectively. 

Figure 4. Dipole power spectra |P(a;)p (arbitrary units) for the excited cluster after 
one coUision (red sohd hne) and after two collisions (black dashed line). Frequency is 
shown in a.u. 
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